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Abstract 

The long-range pt and multiplicity (n) correlations in high-energy nuclear col- 
lisions are studied in the framework of a simple cellular analog of the string fusion 
model. 

Two cases with local and global string fusion is considered. The pt~n and 
n-n correlation functions and correlation coefficients are calculated analytically 
in some asymptotic cases using suggested Gauss approximation. 

It's shown that at large string density the pt~n and n-n correlation coefficients 
are connected and the scaling takes place. The behavior of the correlations at 
small string density is also studied. 

The asymptotic results are compared with results of the numerical calcula- 
tions in the framework of proposed cellular approach. 







1 Introduction. 



The colour strings approach |], 0] is widely applied for the description of the soft part 
of the hadronic and nuclear interactions at high energies. 

In the frame work of this approach the string fusion model was suggested in papers 
Pfl. Later it was developed |§-[H] and applied for the description of the long-range 
multiplicity and pt correlations in relativistic nuclear collisions 

In the paper we have formulated some simple cellular analog of the string fusion 
model, which enables explicit analytical calculations of the correlation functions in some 
asymptotic cases and can simplify calculations in the case of real nuclear collisions. 

In that paper we have checked up the assumptions of the cellular approach and the 
validity of a suggested Gauss approximation in the simplest (no fusion) case when the 
explicit solution of the model can be found. 

In the present paper in the framework of proposed cellular approach we calculate 
Pt-n and n-n correlation functions and correlation coefficients in two cases: with local 
and with global string fusion. 

The calculations are done both numerically and in some asymptotic cases analyt- 
ically using the Gauss approximation. The results of both calculations are in a good 
agreement, which proofs the validity of proposed Gauss approximation. 

At large string density 77 ^> 1 the connection between the pt-n and n-n correlation 
coefficients are found in both cases: with local or global string fusion. At that for 
the correlation coefficient the scaling takes place. It depends only on one combination 
Uo/yfrj of the variables 77 and //q (/xo is the mean multiplicity emitting by a single 
string). 

The paper organized as follows. In the next section we recall the formulation of the 
cellular approach in the case with a local string fusion. 

In the section 3 we develop Gauss approximation at large string density 77 in the 
local fusion case. 

In the section 4 the pt-n and n-n correlation coefficients are calculated at large 77 in 
the case with local string fusion and connection between these two coefficients is found. 

In the section 5 the pt-n and n-n correlation functions and correlation coefficients 
are calculated at large 77 for "homogeneous" situation (a constant mean density of 
strings) in the case with local fusion. It's shown that in this situation the pt-n and 
n-n correlation coefficients become equal and the /io/y^-scaling takes place. The 
obtained results are compared with the results of numerical calculations using formulas 
of section 2. 

In the section 6 we recall the formulation of the cellular approach in the case with 
a global string fusion. The exact closed formulas for the pt-n and n-n correlation 
functions in this case are obtained. 

In the section 7 we develop Gauss approximation at large string density 77 for the 
global fusion case. The pt-n and n-n correlation functions and coefficients for this case 
are calculated and compared with the results of numerical calculations using formulas 
of section 6. 

The behavior of the correlations at small string density is studied in Appendix A. 
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2 Cellular approach to the local string fusion. 



Let us recall the formulas obtained for this case in |Tl| . We consider the collision of 
nuclei in two stage scenario when at first stage the colour strings are formed, and at the 
second stage these strings, or some other (higher colour) strings formed due to fusion 
of primary strings, are decaying, emitting observed particles. 

In principle, one can consider two types of fusion. The case with a local fusion 
corresponds to the model, where colour fields are summing up only locally and the 
global fusion case corresponds to the model, where colour fields are summing up globally 
- all over the cluster area - into one average colour field, the last case corresponds to 



the summing of the sources colour charges. (In section 5 of [11] we have referred to 
these cases as A) and B) correspondingly.) 

In the transverse plane depending on the impact parameter b we have some inter- 
action area S(b). Let us split this area on the cells of order of the transverse string 
size. Then we'll have M = S(b)/a cells, where o"o = tttq is the transverse area of the 
string and r$ Pd 0.2/m is the string radius. 

In the case with a local fusion the assumption of the model is that if the number of 
strings belonging to the z-th cell is rji, then they form higher colour string, which emits 
in average fioy/Vi particles with mean p\ equal to p 2 ^/rfii compared with fi particles 
with (pi) = p 2 emitting by a single string. 

If we denote by rij and - the number and the average number of particles emitted 
by the higher string from i-th cell in a given rapidity interval, then 

m = noVvi (i) 

From event to event the number of strings rji in i-th cell will fluctuate around some 
average value - rj i . Clear that in the case of real nuclear collisions these average values 
Tj i will be different for different cells. They will depend on the position (s) of the i- 
th cell in the interaction area (s is two dimensional vector in transverse plane). To 
get the physical answer we have to sum the contributions from different cells, which 
corresponds to integration over s in transverse plane. 

The average local density of primary strings r\ i in the point s of transverse plane 
is uniquely determined by the distributions of nuclear densities and the value of the 
impact parameter - b. They can be calculated, for example, in Glauber approximation. 
We'll do this later in a separate paper. In present paper we consider that all r\ i are 
already fixed from these considerations at given value of the impact parameter - b. 

Let us introduce some quantities, which will play important role in our considera- 
tion: 

M M 

# = 5>, n = y,% 

1=1 8=1 

M M 

r = Ev^. r = Y J \fii (2) 

i=i i=i 

then clear that N is the number of strings in the given event and iV is the mean number 
of strings for this type of events (at the fixed impact parameter b). 

To go to long-range rapidity correlations we have to consider two rapidity windows 
F (forward) and B (backward). Each event corresponds to a certain configuration 
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{r] 1 , ...,r] M } of strings and certain numbers of charged particles {n 1 ,...,n M } emitted 
by these strings in the forward rapidity window. Then the total number of particles 
produced in the forward rapidity window will be equal to rip: 



M 



n F = J2ni (3) 



i=l 

The probability to detect rip particles in the forward rapidity window for a given 
configuration {r^, ...,i] M } of strings is equal to 

M 

P{ Vv ...,V M }(n F ) = E tnr^TnHPntfa) ( 4 ) 

f»lr-,« M ) i=1 

where p (rij) is the probability of the emission of n, particles by the string rji in the 
forward rapidity window. By our assumption ([!]) 

oo 

rk = E ^PrjM^ = VoVVi (5) 

rtj=0 

If we denote else by W(rji, ...,i] M ) the probability of realization of the string con- 
figuration {r)i, ...,r] M } in the given event, then the average value of some quantity O 
under condition of the production of np particles in the forward window will be equal 
to 

,V M }' n F 

[ K s hl , m1 w(%,,%)\ , m1 (%) 1 } 

One has to omit in this M-fold sums one term, when all rji = 0, which corresponds 
to the absence of inelastic interaction between the nucleons of the colliding nuclei (see 
details in Appendix A). 

If the O in the number of particles produced in the backward rapidity window n B 
in the given event, then we have to use for (n B ) n correlations: 

At 

( n fl){ii,....i M W =/ f oEv / ^=W (7) 

i=\ 

If the O in the mean squared transverse momentum of particles produced in the back- 
ward rapidity window p\ B in the given event, then we have to use for {p^ B ) n correla- 
tions: 

(p 2 tB) {Vl ,..., VM },n F = E t m ^ P'VVI = P 2 m 1 1 = P 2 - (8) 
i=iLi=lV'/i ^1=1 v'/« ' 

Later we'll assume that numbers of primary strings in each cell rji fluctuate inde- 
pendently around some average quantities f] i uniquely determined by the distributions 
of nuclear densities and the value of the impact parameter - b (see above), then 

M M 

W(t]i,-,Vm) = I[w(Vi), J2vMVi)=Vi (9) 
i=i i=i 

For clearness we'll sometimes address to a simple "homogeneous" case, when all 
rj i (but not the rji, which fluctuate from event to event!) is equal each other in the 
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interaction area l] i = 77 (a constant mean density of strings). The parameter 77 coincides 
in this case with the parameter 77 used in the papers || || |§ and has the meaning of 
the mean number of strings per area of one string (rj = (mean string density) x cr ) . In 
general case the parameters r\ i have the same meaning, but with mean string density 
depending on the point s in the transverse interaction plane 
(7^ = (mean string density in the point s) x cr ). 

As it was shown in |T]J if we assume else the Poissonian form of p^iyii) (p a (x) is 



the Poisson distribution with x = a): 



P^-P^in^e-^J^T (1Q) 

then we find the Poissonian distribution for 

P {Vl ,..., nM }(n F ) = P^ oEiVm (n F ) (11) 

with (n F )^ j )?? -j. = /ioT,i y/Wi = Por = (n F ) r and o 2 rF = (n F ) r = /x r, which we can 
replace at large /x J2i y/Vi by Gauss distribution: 



(n F -(n p ) r ) 



p <- =^r (12) 



It was also shown in the section 6 of fill that if we assume the Binomial form of 



p Vi {rii) then we find the Binomial distribution for ^ y(n F ) with (n F )^ v ^ y = 

PoJ2i yffji = Po r = { n F)r an d a rF = ( n F)r(^~^) = A i o r (l — A), where A — > corresponds 
to the Poisson limit and A — > 1 corresponds to the case when each higher string emits 
fixed number of particles - rii = rii = fioy/vl in each event. 

Moreover appealing to the central limit theorem of the probability theory one can 
state that at large M we'll have formula ( |12|) for any type of p (rii). 

3 Gauss approximation at large string density for 
the local fusion. 

Now we'll evaluate J2{r tl ....,ri M } & t the condition that all rj^ ^> 1. At this condition we 
can use Gauss approximation for each w(r)i) 

w( Vi ) = - 7 =^e (13) 

V ZTT<T„. 



with o~. = 7^(1 — X v ), where again — > corresponds to the Poisson limit and A^ — > 1 
corresponds to the case with a fixed number of strings N = N = J2fli Vi ( see the 
section 6 of 0]). 



As in [111 we have then 



f dr) x ...dr] M r^e 



-<P(.Vi,n F ) 



J dr] 1 ...dr) M ^e Wfc.»*0 
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and 

Recall that JV = r)i and r = YllU VVi- Here 

fri 2^(1 - A^) 2/i r(l - A) 

Further we take out factors before exponent at the point, where <p is minimal, after 
that the rest integrals in the numerator and in the denominator are being reduced, and 
we find 

( n B)n F = V0 r * (17) 

and 

N* 

(p"b)u f =P 2 — (18) 

The N* and r* are the values of N and r in the point {77*, ■■■,'rf M }, where ip(r]i,n F ) is 
minimal: 



d<p(rn,n F ) 
drji 

This leads to the system of equations: 



(19) 



Recall that r* = y/rj? and 

" = ^ < 21 > 

For the meaning of k see the end of the previous section and the section 6 of flllf , for 
both Poissonian distributions k — 1 and k is the relative width of the p(rij) and ty(7/i) 
distributions in other cases. The fl2"0| ) defines rj* as function of n F . 
Introducing short notations 

/ = ^ = 7% a l = ^ (22) 
/i r (n F ) 4^ 

we can rewrite (EDI) as 



^ 3 " ^ = «i (/ 2 ^ - l) ( 23 ) 

where r = E^i v^i> r * = E^i^V^ and ^* = ^i z ili- The (H defines * as 
function of /: ^ = ;%(/). Then we can calculate (n B ) n and (Pt#) n as a function of 
n F using ( |TT[ ) and (|TE|). 
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4 Pt—n and n—n correlation coefficients at a large 
string density for the local fusion. 



The correlation coefficients are defined in the same way as in the section 4 of pi 



(24) 



dn F '"f=K) 
and 

or for "relative" quantities 



_ \ Bl ripl \ Bl . . . 

- dn F /(n F ) ln F=K) 1 u; 

and 

^_ d(rf B ) nF /(p* B ) 

1 dn F /(n F ) K=(n F ) y > 

(Note the same definition of pt-n correlation coefficient /3 in [|12] (see formula (44) in 
p2||), see also remark in Appendix B). 

In short notation using fll"?!) and (0) we have: 

1 dr* 

»=FW^ < 28 > 

and 

T_ d{N*/r*) 

P ~N df '' =1 [J) 

We can't solve the equations ( p3|) for to find z,-, = Zi(f) explicitly, but to calculate 
the correlation coefficients we need to know only z'^l) = dZi ^p |j =1 , which can be done 
explicitly. 



We see that at / = 1 (|23|) has the obvious solution: 

f = l, Zi = l, r* = r, N* = N (30) 



We need to calculate z[(f) only at / = 1. Differentiating (|2"3] ) on / and using then 
again ([30]) we find 

m = a _f (31) 
4r + HqhM 

with di = fi n/(4:y/^). Then 

t Idr* 1 *C , rr liqk 

b = --TT \ f =i = -2^ z i 1 )yVi = — j_ A -, A/r 32 
r df 3 r ~{ v yU K + 4r/M 

and 

1 dN* _ldr* ldN* _r 

Jf df r df " TV df ° 



Using 

dN' " , 2^r* 

lT l '- = 2 S 2 ' (lk = ^T* (34) 

we have 

\NM J hqk + Ar /M \NM ) v ; 

We see the connection between pt-n and n-n correlation coefficients. Note that due 
to obvious inequality: 

/ m \ 2 M 

EV^ < M J2% (36) 

\i=l J i=l 

we have r 2 < MiV and hence always (3 < b. 

Clear that at equal rj i = 77 we have f = My^rj, N = Mn, r 2 = NM and 

]3 = 6 = ^ - (37) 

From ([32]) we see that n-n correlation coefficient is always positive. Can pt-n 
correlation coefficient be negative? Let us consider nonhomogeneous situation when 
jj. = n + at i = 1, Mi, % = r}_ at i = M 1 + 1, M, M 1 ~ M and 77+ > ?7_ > 1. 
Then r = M lv ftjT + (M - Mi) v /7tZ « M ly /Tf£, N = Min + + (M — Mi) 77- « M1774. 
and we have 

? « (2^± - l) 6 (38) 



M 

We see that at M x < M/2 we can have ~j3 < 0. 

5 The fio/rj 1 / 2 -scaling at large string density. 

Let us consider for clearness homogeneous case, when all r\ i is equal each other in the 
interaction area 7} i — 77 (a constant mean density of strings). In this case we can 
explicitly calculate at large string density rj not only the pt-n and n-n correlation 
coefficients, but also the corresponding correlation functions for the version with a 
local string fusion. 

We have seen in the end of previous section that in this case the coefficients for n-n 
and pt-n correlations defined as (p6|) and ([27]) are coincide: 

? = & = — ^ = "^7 (39) 
Hon + a + 1 

where a = /i /t/(4y / r/). 

In this homogenous case Tj i = 77 we can also calculate the correlation functions 
(PtB)n F an d ( n B)n F a ^ an y n F- Due to symmetry, the system of equations (^3|) have 
symmetrical solution Zi — z and can be reduced to one equation 

z 3 -z = a(^-l) (40) 
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because r = M^rj, r* = zM^Jf) and N* = z 2 Mi] with 



/ioM^r/ (n F )' 4^7/ 1 ; 

The (1^) defines the function z = z(f) and then using ([17]) and (|l"8"D we can calculate 

(n B )n F = = ^M^z(f) = (n B )z{f) (42) 

and 

N* 

(pIbK = V 2 — = pVW) = ( P 2 tB )z(f) (43) 
So we have at any n F = (n F )f: 



{n>B)n„ (PtB) 



= z(f) (44) 



(n B ) {p 2 b) 

From (|39|) and ( |40|) we see that in this homogenous case at large string density rj 
there is an remarkable scaling. The pt-n and n-n correlation coefficients and correlation 
functions depend only on one combination a = ixqk/ {A^/ff) of parameters. 

We present the correlation function z(f) (|4*4| ) in Figs. and the correlation co- 
efficient P — b ( |5SD as function of 77 in Figs. |]-| (the solid lines). We present also in 
Figs. the results of our direct numerical MC calculations of the pt-n and n-n cor- 
relation coefficients in the local fusion case based on formulas (§-§]) (empty and filled 
points correspondingly). 

We see that in the case with local fusion at small string density we have large n-n 
correlations (the same as in the case without string fusion [|lTJ) and no Pt~n correlations. 



(The analysis at very small values of rj < 1/M see in Appendix A.) 

At large string density in the homogeneous case the pt-n and n-n correlation coef- 
ficients become equal and the /1q / y^-scaling takes place. We see also that in this limit 
our Gauss asymptotic is in a good agreement with results of the numerical calculations 
and M-independence takes place. 



6 The global fusion at large string density. Exact 
solution. 

In this case at first stage we also have M = S(b)/ao cells (like in the case with local 
fusion) with rji, i = 1, M fluctuated around T] i . Then (unlike the local fusion case) we 
have to find average r/ c = J2i Vi = Ji f° r given event, and then to generate particles 



from one cluster with average multiplicity equal to u cv /7^ = HoM^/rfc = fi Q M^jN/M 
= /iqa/ MN. The general formulae for this case was obtained in the section 5 of [|TT|j: 



E {??1 ,.., %f} ^(r /l ,...,r /M ) Vv/¥r _(n F ) 

where \i c = [IqM with M = S(b)/a . The assumption r\ i » 1 is essential, as only in 
this situation we can consider that the transverse area of the cluster AS is equal to all 
interaction area S(b) {b- impact parameter). 
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The (O)^ v } n is the rates of the backward production from configuration 
{t7 1; ...,r] M }. We have to use for (n B ) n correlations: 



and for (p\ B ) n correlations: 



\ 



1 EVi (46) 



M 



I 2 \ _ 2 / — _ 2 



\ 



We see that the difference with the case of local fusion consists in replace v 7 ^ 



As a consequence calculations in the case of global string fusion are much 
more simple, as we can reduce all sums J2{ Vl ,...,ri } to one sum J2ni as m the no fusion 
case (see section 3 in |[1|). So in the global fusion case we can write simple formulas: 

{ " bK = e»«wp m w*> (48) 

and 



Y.„W(N)p K ^{n F ) 
where W(N) is given by the formula 

W(N)= J2 (50) 

We see that in the case of global fusion (one cluster at large TJ i with area AS being 
equal to all interaction area 5(6)) n-n and pt-n correlations are connected 

2 P 2 I \ ( n B)n F (PtB)n F , . 

<»->., = ^<» B >., - "kT = w (51) 

Note that unlike the local fusion case in this case we find this result without any 
assumptions on the properties of p(n F ) and w(j]i) distributions and for arbitrary (even 
nonequal) r) i . 

Clear that in this case the results can depend only on mean number of strings N 
and on combination hm'- 

N = Y,% m = »aVM (52) 



Below we'll calculate numerically the correlation functions on formulas ( B8[ ) and 
( fi9|) , but at first we would like to find explicit formulas for global fusion case in Gauss 
approximation. We'll see that results really depend only on one combination of the 
variables (p2"D , namely on 

PM (53) 



and the scaling takes place as in the case with local fusion. 
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7 Gauss approximation for global fusion at large 
string density. 

Acting as in the no fusion case in section 4 of (see also calculations in the Gauss 
approximation in the section 3 of the present paper) we find 



(n B )n F = /iA/ViV* = fioVMVN* (54) 

or keeping in mind (CT) 



( n B)n p _ {PIb)u f jN*_ 

(n B ) ~ (P? B ) "ViV 
The N* is the value of N at which the function 



(55) 



^ F) 2N(1-X n ) 2(jl m VN(1-\) 



gains its minimum. In short notations 



N* n F n F [x m k h k , . 

z = W^, f = = = - — -, a— — = = — . (57) 

V n MM Viv M a^/n aJn/m 



we find the equation 



z 3 - z = a I J — - 1 1 (5* 
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2 



-2 



which defines the function z = z(f) and then using (|55|) we can calculate correlation 
functions 



( n B)n„ (PtB)n 



n 



Mo = *(/) (59) 



and correlation coefficients for the global fusion case 

? = W(1) = — = ^ ^ = (60) 

We see again that in Gauss approximation there is the same remarkable scaling. The 



MOW 



pt-n and n-n correlations depend only on one combination of parameters: a — ,= — . 

iy/N/M 

Note that unlike the local fusion case in this case we find this result for arbitrary (even 
not equal) fj^ 

In the homogeneous situation all r\ i = rj and we have 

N = Mn, a = ^- (61) 
"4^ 

and 

3 = 5=^-=^- (62) 
/io« + 4^ a + 1 

We see that in the homogeneous situation in Gauss approximation the results with 
local and global fusion coincide (as we have expected in section 5 of |i~l| ). We have the 
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same equations ( f40|) and (|58| ) with the same value of parameter a ( f4"T|) and (|61|). Note 
that in the no fusion case we have had very different equation for z(f) (see section 4 

of H) ; 

Unlike the local fusion case in the global fusion case we can control the validity of 
Gauss approximation making calculations on exact formulas ((4*8|) and (f49|) at different 
values of M. 

Along with the correlation function z(f) (fPj), ( |59D in Figs, ffl^ and the correlation 
coefficient (3 = b (p9|),(|60D in Figs. [3|-[5], calculated on our scaling formulas, which in 
Gauss approximation are the same for local and global fusion (the solid lines), we 
present at the same pictures the results of exact calculations in the global fusion case 
on formulas (|48|) and ( f4~9D at different values of M (the dotted and dashed lines). 



We present also in Figs. [||| the results of our direct numerical MC calculations of 
the pt-n and n-n correlation coefficients in the global fusion case based on formulas 
(half-filled squares). 

We see that in this case the Gauss approximation works very well and the fio/y/rj- 
scaling is not an artifact of this approximation. More over along with the /io/ y^-scaling 
at large t] we have also M-independence for correlation coefficients (3 and b (see Figs. 
[5]) starting very early (from M = 4). 

8 Conclusions. 

In conclusion let us compare the results obtained in the present paper in the case with 



string fusion with results obtained in JTT[ in the case without string fusion. 



We have obtained in the paper |l I]] in the case without string fusion: 

1. for n-n correlations: b = -2- with a = iinK 

2. for pt-n correlations: (3 = 

In the present paper in the case with the global string fusion and in the local 
fusion case for homogeneous situation (r/j = if) we find at large rj: 



1. for n-n correlations: b = -tt with 

a+1 



2. for pt-n correlations: (3 = b = with the same a 

We see that with fusion the n-n correlations became weaker, but now as compensation 
we have the pt-n correlations of the same strength. We see also fio/y/rj- scaling in this 
case. 

For nonhomogeneous situation (different rjf) in the case with local string fusion 
we have find at large r\{. 

1. for n-n correlations: b = with a = fj, K/(Ar/M) 

2. for pt-n correlations: (3 = — ljb and hence (3 <b 
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M M 



with N = Vi an d r = ^ \/v 



i=l i=l 



As we have demonstrated above (see ( |33"D and (pq)), this leads to /3 smaller than 6: 
(3 <b. It's possible situation (j38l), in which /3 < 0. 

At small string density as it's shown in Appendix A the two types of limit at n — » 
can be studied. 

1. If one keeps M = const, then we have iV — > 1 (because the configurations with 
iV = 0, are not considered as events) and we have nor pt-n nor n-n correlation. 

2. If one keeps N = const, then Mi] = const and M — > oo, hence the strings will be 
far separated in transverse plane and we'll have the same results as in the case 



without string fusion (TT| 



Note that the results obtained in our cellular approach are in a good agreement 
with the results obtained in the framework of the real string fusion model taking into 
account detail geometry of strings overlapping |13 |. 
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Appendixes: 

A Correlations at small string density. 

In this appendix we calculate the correlation functions and the correlation coefficients 
at small string density in the case with local string fusion. (In the case with global 
fusion this limit has no physical sense (see discussion in the section 5 of [|TTlD ). 

For clearness we consider " homogeneous" case, when all Tj i is equal each other in the 
interaction area r/j = n. Then in each cell i (i = 1, M) the r]i fluctuate around this 
mean value according to Poisson law {p a (x) is the Poisson distribution with x = a): 

w{m) = P v (Vi) = (63) 
We'll assume also the Poissonian form of p Vi (ni): 

= W^) - e-^i^r (64) 

then we have the Poissonian distribution for 

P {Vl ,..., VM }(n F ) = P, oEiVm ;(n F ) (65) 
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with (n F ) {r)ii Vm] = ^oEi ^Wi = Vor = (n F ) r . 
Then we find from (IBT) for n-n correlations: 



^oS / {?7l „„,, M} r(n^fa))p w? .M 



and for pt-n correlations: 



P 2 Ef„ , w „i ~ (niW(?7i))p„ or (n F ) 



Recall that 



tBnp E' {ril ,..., VM} (UMVi))p, or (n F ) [ ) 

M M 

JV = ^ 77; and r = Vvl 

i=l i=l 

The JV is the number of strings in the given event. 

One has to omit in all M-fold sums in (|66f) and (|67| ) one term, when all r]i = 0, 



which corresponds to the absence of inelastic interactions between the nucleons of the 
colliding nuclei (we denote this fact by J2')- 

The probability P(n F ) to detect n F particles in the forward rapidity window, which 
enters denominators of ( |66|) and (67), is 



P(n F )=C J2 ' (IW^)J P» rM 



(6£ 



where from normalization condition we have 



° l-w M {0) l-e- M " (69) 

Clear that this factor C is canceling in the numerator and in the denominator of 
(|B5|) and (p7|), but if we calculate the mean number of strings JV, we find 

and for the (n F ) at small 7/ 1 we have 

M = J2 n F P M = C^OF E' [H w (jh) ) r = NF , ^-Mt, = ^OF^ ( 71 ) 

Because for any w > we have J2 m r ]i w (' r )i) = V + 0{rf) at 77 — >• 0, as the main 
contribution comes from the term r/j = 1. 
There are two possibilities when 77 — ^ 0: 

1. M = const and then Mr/ -> and JV ->• 1 (see fl70|)) 

2. JV = const and then Mr] = const ((see (fT0|)) and M — > 00 
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We'll investigate both these possibilities. 

The first one means that in the limit we have N = N = 1 (because the configura- 
tions with N = 0, are not considered as events). Clear that in this situation we'll have 
nor pt~n nor n-n correlation, as we'll have no fluctuations in the number of strings (see 
discussion in the end of the section 4 of [0). Detail calculations are presented below. 

In the second case with N = const we'll have the fluctuations in the number of 
strings N, but in the limit 77 — > we'll have M — > 00, the strings will be far separated 
in transverse plane and the strings fusion will plays no role. So in this case we'll have 



the same results as in the no fusion case, considered in fill : large n-n correlation 
with a correlation coefficient equal to b = /x / (/x + 1) and no Pt~n correlation. (See 
calculations below.) 

Detail calculations. Let us evaluate the M-fold sums in (BB) and ( |67j ) at n — > 

keeping all terms of order (Mn) k , (Mt]) k r] with any k and omitting all terms of order 
(Mn) k r] 2 and higher. 

The terms of order (Mrj) k originates from the summands in (|66|) and (67) with 
T] il — ... — r)i k — 1 and other rji = 0. The terms of order (Mn) k n originates from the 
summands in (|66D and (|67j) with 77^ = 2, r) i2 — ... — r\i k — 1 and other rji = 0. Keeping 
this into mind we find for P(np): 

P(n F ) = C(G + Id) (72) 



and for n-n correlations: 



N + gjVi 

\ n B)n F = VoB Go + RGi ( 73 ) 



<&>.. = P 2 ^M- = P 2 fl + (74) 



For pt-n correlations we have: 

l ^ = " Go + f Gi * K ' 2 G, 

Af 

G = Po = EcWp^(M (75) 



Here 



fc=i 

M 

fc=i 
fe=i 



,1/ 



Mo(fc+7)^J 



fc=l 
M 



fe + 1 



where /i = /i 0F and 7 = — 1 . 

Note that at M ^> 1 and Mr/ = const we have for P(uf) 



P(n F ) = Ce - M ^(G + V^7 Gl) (76) 
nJ. 2 
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and for n-n correlations: 



N + 2g-Mo7 j /\T 1 



G + U-wd 



(77) 



For pi-n correlations we have: 



(Pw); 



P ^ 



Here 



Go = Pq 

k=l 



Gc 



fc=l ft/ - 

x Jk 
n F + 1 — 
k\ 



(78) 



(79) 



k=l 

oo 

iVi = E( fc + 7) n - +1 



d k 



k=l 



(k-i)\ 



(k-iy. 

where d = Mrje"^, 

For the control of calculations we have used also the following explicit formulas at 
n F = 0,1,2: _ _ 

G (0) = P o (0) = e d - 1 (80) 

iVo(O) = G!(0) = G (l) = P (l) = de d 
N 1 (0) = 'G 1 (1) = (d + V2)de d 

iVo(l) = Go(2) = P (2) = (rf+l)cie d 
iVi(l) = Gi(2) = (d 2 + d(l + 2^2) + 2)de d 
= (rf + 2)rfe d 
JV (2) = (d 2 + 3d+l)cfe d 
]Vi(2) = (d 3 + d 2 3(l + V2) + d(7 + 3V2) + 2y/2)de d 
Pi (2) = (d 2 + d(3 + V2) + 2V2)de d 
The results of the calculations using the asymptotic formulas (|73| -[75|) in the first case 



(M = const, 7] — > 0, Mr? — > and iV — > 1) are shown in the Figs. |6|-(7] (lines) together 



with results of direct MC numerical calculations using formulas (|66|) and (|67|) (points). 
As we have expected in this case both n-n and pt-n correlation coefficients go to 
zero when N ~ 1, i.e. at n < 1/M. Remember that we have M-independence for 
the correlation coefficients at large r\. Now we see that in this limit it disappears at 
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n < 1/M. This is also the reason for nonlinear dependence of the correlation coefficients 
on T) in this region which one can see in Figs. [p-|7] for /^o = 4. 

We see also that using asymptotic formulas (|?3|-|T5|) we can calculate n-n correlation 



coefficients in wider region of small 77, than pt-n correlation coefficients, because we 
have the contributions of order (Mn) k and (Mn) k n for n-n correlations and only first 
non-trivial contribution of order (Mn) k n for pt-n correlations. Note the very good 
agreement between results of the calculations on the asymptotic formulas (|73|-[75D (lines 
in the Figs. ^-0) and the results of direct MC numerical calculations on formulas (|66|) 
and (|67|) (points in the Figs. §-0). 

In the second case when n — > we keep N = const and then due to flTOD Mn = const 
so M — > oo. We can use formulas (ff7| - f79|) with d = const, then in the limit n — > we 
find for n-n correlations: 

/ v _ SgLi kp Mv (k) PlMQk (n F ) 

{n B ) np - ^ - »ob E oo =i p^p^) ^ 

Here we multiply both the numerator and the denominator by c~ Mr] [i^ F jn F \. 

Recall that p a (x) is the Poisson distribution with x — a, then we see that formula 



(pl|) coincides with the formula for n-n correlations obtained in the paper |ITJ in the 
case without string fusion (see formula (23) in JTTJ). So we'll have the same result for 
the n-n correlation coefficient b = /U /(/i + 1) as in the no fusion case (see Figs. ^-0). 
For pt-n correlation coefficients in this limit we find from ( |7H|) and (|79| ) with d = 

const: 

(p 2 tB ) n =p 2 (l + 0( V )) (82) 



So we have no p t -n correlation as in the case without string fusion [If . 



B On the difference between (n B ) and (n B ) n = ^ n y 



It's possible instead of (|26| ) and (|27|) to use the following definitions for the correlation 
coefficients: 



and 

where 
and 



dn F j (n F ) 

-_ d(p 2 B)n F /(PtB)(n F ) , ( . 

P~ dn F /(n F ) '^=<n,> 

( n B)(n F ) = ( n B)n F = (n F ) (85) 
(Pt B ){n F ) = (Pts)n F = (n ) (86) 



Clear that 

l B) = Y,P( n F)( n B)n F ~ ( n B)n F = (n F )J2 P ( n F) = ( n B )n F =(n F ) = M (n F ) (87) 



n, 



l F~ V'FI 

1p Tip 



and the same for (p 2 B ) {nF) = (p 2 B )n F =(n F ) ~ (Pib)- 

In our Gauss approximation these two types of quantities coincide with each other. 
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with fusion, muO/sqrt(eta)-scaling 
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Figure 1: The pt-n and n-n correlation functions. Solid line - the Gauss approximation 
for local and global fusion at the value of scaling variable /in/y^ = 1- Dashed and 
dotted lines - the results of exact calculations in the global fusion case at different 
values of /i , V an d M. The Ho/ y/rj- scaling and M-independence. 
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global, mu0=4,eta=2,M=4 
global, muO=4,eta=2,M=128 
local and global, mu0/sqrt(eta)=2.83, Gauss 
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Figure 2: The same as in Fig. [T], but at different values of /xq and r\. 
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correl.coeff., with string fusion, mu0=1 



local & global, Gauss approx.,[muO/(muO+4sqrt(eta))] 

global, M=1 28 
global, M=4 
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global M=30 n-n & pt-n 
local M=30 n-n 
local M=30 pt-n 
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Figure 3: The pt~n and n-n correlation coefficients at /in = 1- Solid line - the Gauss 
approximation for local and global fusion b = (3 = Uq/(uq + ^\/v)- Dashed and dotted 
lines - the results of exact calculations on the formulas fl48[ - fl9| ) in the global fusion case 
at different values of M (half-filled squares - the same by means of direct numerical 
MC calculations on formulas (p5| - [47|) ) . The empty and filled points - the results of 
direct numerical MC calculations in the local fusion case based on formulas (0-^) for 
the pt-n and n-n correlations correspondingly. The filled circle - the n-n correlation 
coefficient b = Uq/{uq + 1) in the case without string fusion [ITT . 



correl.coeff., with string fusion, mu0=2 
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local & global, Gauss approx.,[mu0/(mu0+4sqrt(eta))] 
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global, M=4 
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local M=30 n-n 
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Figure 4: The same as in Fig. |], but at uq = 2. 
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Figure 5: The same as in Fig. but at /in = 4. 



n-n correl.coeff. at SMALL eta (with local string fusion) 
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Figure 6: The n-n correlation coefficient at small values of n for /j = 1 and /in = 4. 
The lines - results of the calculations using the asymptotic formulas (|73|) and (|75| ) at 
M = const (N — ► 1). The points - results of direct MC numerical calculations using for- 
mula fl66l). The arrows show the value of the n-n correlation coefficient b = /in/(/Jn + l) 
in the case without string fusion 



11], which corresponds to the limit N = const 
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pt-n correl.coeff. at SMALL eta (with local string fusion) M=30 
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Figure 7: The pt~n correlation coefficients at small values of n for /x = 1 and /io = 4. 
The lines - results of the calculations using the asymptotic formulas ([74]) and (|75| ) at 
M = const (N — > 1). The points - results of direct MC numerical calculations using 
formula 
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